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f— ' . Abstract 

We investigate the effects of pair creation on the internal geometry of a 

£> ■ black hole, which forms during the gravitational collapse of a charged massless 

q^ | scalar field. Classically, strong central Schwarzschild-like singularity forms, and 

^ir \ a null, weak, mass-inflation singularity arises along the Cauchy horizon, in such 

a collapse. We consider here the discharge, due to pair creation, below the event 

horizon and its influence on the dynamical formation of the Cauchy horizon. 

^ { Within the framework of a simple model we are able to trace numerically the 
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collapse. We find that a part of the Cauchy horizon is replaced by the strong 
space-like central singularity. This fraction depends on the value of the critical 
electric field, E CT , for the pair creation. 

1 Introduction 

The well known exact solution of the coupled Maxwell-Einstein equations outside the 

spherically-symmetric matter distribution is the Reissner-Nordstr0m solution. The 

analytic extension of the Reissner-Nordstr0m metric has rather exotic properties. The 

black hole's interior contains Cauchy horizons, time-like singularities and tunnels to 

other asymptotically flat regions. 

Recently, it has been shown both in perturbative analysis and by solving the full 

non-linear problem that a Cauchy horizon inside a charged black hole is transformed 

into a null, weak singularity |TJ, g, |3], [|, |5], 0], The Cauchy horizon singularity is 
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Figure 1: Penrose diagrams, (a) and (b) are, respectively, the Schwarzschild and the 
Reissner-Nordstr0m eternal (pre-existing) spacetimes , (c) and (d) are ones expected to 
form during a dynamical collapse of a neutral and a charged matter, respectively. The 
actual, dynamically calculated, space-time geometry for a charged gravitational collapse, 
(e) is dramatically different from that depicted on (d) and it resembles (in broad terms) the 
Schwarzschild space-time. In all the figures thick solid lines represent central singularities 
(space-like for the Schwarzschild manifold or time- like for the Reissner-Nordstr0m one). 
The weak, null singularity in figure (e) is represented as a thick dashed line, while thin 
dashed lines describe various horizons. 



weak in the sense that an infalling observer crossing it experiences only a finite tidal 
deformation |3|, |J. However, the curvature scalars diverge along the Cauchy horizon, 
leading to an unbound growth of the internal mass-parameter, a phenomena known 
as the mass-inflation ||. The earlier studies were done on the pre-existing (eternal) 
Reissner-Nordstr0m space-time. Hod and Piran J7|, have demonstrated explicitly that 
mass-inflation takes place also during a dynamical charged gravitational collapse. The 
dynamical space-time is drastically different from the analytically extended Reissner- 
Nordstr0m manifold: it resembles more the Schwarzschild one. The Penrose diagrams 
of the various space-times are depicted on Figure |l|. 

This is so far the classical picture. Our goal, here, is to consider quantum effects 
and to investigate the influence of pair creation in strong electric fields on charged 
gravitational collapse. Specifically, we are interested in the effects of pair creation 
on the inner structure of the black hole that forms in such a collapse. This goal 
differs from previous works on the subject. Pair creation was mainly considered in 



the external region of a pre-existing black hole's space-time, outside the event horizon 
H |TD], [II], |T2|. It has been shown that the produced particles rapidly diminish the 
charge of a black hole as seen by an external observer. 

Particles creation takes place, however, also in the inner region of a black hole. 



Novikov & Starobinskii |L3[ and Herman & Hiscock [JJj studied the inner geometry 
and the stability of a Cauchy horizon in the pre-existing Reissner-Nordstr0m space- 
time influenced by the pair creation effect. The model in Ref. HM assumes the 



instantaneous disappearance of the electric field, when the pair creation takes place 
inside the event horizon along r cr = const hypersurface. The Reissner-Nordstr0m 
patch of the space-time exterior to this hypersurface is glued along the r cr = const 
hypersurface to an interior Schwarzschild patch. In this model the Cauchy horizon 
does not exist. The model in Ref. fl3| assumes an initially Reissner-Nordstr0m 
geometry and allows the evolution of the electric field through the back-reaction of 
the created pairs. In this model the initial Reissner-Nordstr0m geometry evolves to 
an uncharged, Schwarzschild-like one - the Cauchy horizon is shown to be unstable 
with respect to the process of pair creation. It should be emphasized that in this 
model the Cauchy horizon is assumed to be stationary tch = const hypersurface. 
This is in contrast to the recent investigations Jl|, ||, ||, |[ |!| |6], |7|], where the Cauchy 
horizon was shown to be non static, contracting null hypersurface. 

Although, the particle production in the intensive electric field is the simplest 
quantum process, its influence on the inner structure of black holes was not studied 
yet in an evolutionary context. The effect of pair creation in the intensive electric 
field is probably most important, when dealing with a formation of the inner structure 
of a charged black holes. This depends, of course, on the parameters of the formed 
black hole. In this work we take a different point of view from |TB|, [TJ]] and explore the 



dynamical picture i.e., we replace the question is the Cauchy horizon stable, with the 
question does it form at all. To address this question, one should consider the collapse 
of a charged self-gravitating matter. The electron-positron pairs are produced in the 
electric field of the collapsing matter. To treat consistently the problem of such a 
collapse, one should take into account the back-reaction of the produced pairs on 
the source's electric field. This is achieved by adding the electric current due to the 
produced pairs as a source to the Maxwell equations: 

***., = 47T J£ ee + 47T J p * airs . (1) 

The charge conservation equation in the situation with pair creation is modified by 



adding the charge source: 

J a . a = T(F a(S ) . (2) 

The stress-energy of the electric current of the produced particles arises from the 
stress-energy of the electric field. The latter is the source of the pairs, by means of 
the energy-momentum conservation. 

To formulate the problem properly we need a back-reaction formalism that in- 
cludes the back reaction of the pairs on the stress energy tensor of the electric field 
that create them. Without this the problem would not be self consistent. Such a 
formalism is not available. Instead we consider here a toy model that utilizes the 
main physical properties of the system - the fact that the pairs limit the electric 
field to a critical value, E CT . We describe this effect of pair creation by introducing 
a nonlinear dielectric constant which prevents the electric field from exceeding E cr , 
the critical pair creating field. In doing so we have ignored the electric current of the 
pairs and their stress-energy We also disregard the contribution from vacuum polar- 
ization, which becomes significant only for the exponentially large fields ||. In spite 
of this simplifying assumptions we believe that this model captures the characteristic 
behavior of the real system. 

In section |2| we consider discharge in a classical space-time and give the motivation 
of the investigation of the influence of the pair creation on the inner structure of 
charged black holes. In section [^ we present the underling physical model. We 
develop the formalism and discuss the applicability of the model. Section [| describes 
our numerical scheme. The results are presented in section |5|. We compare a classical 
charged gravitational dynamical collapse with a collapse with the discharge. We 
summarize our conclusions in section |6|. We use units in which c = G = h = 1. 

2 Discharge in a Classical Charged Space-time 

When pairs are created in an asymptotically flat region one of the particles, having 
the same charge as the field's source, is repulsed from the body and escape to infinity. 
The another member of the pair is attracted to the body, decreasing its charge. This 
occurs, for example, to pairs created in the field of a charged black hole, outside 
its event horizon. The black hole discharges rather quickly, until its external field 
becomes subcritical || [K], [0], |T2|1 . 



A more interesting situation occurs when a significant pair creation takes place 
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Figure 2: Pair creation inside the event horizon. Here we depict the Penrose diagram of 
the Reissner-Nordstr0m space-time. Pairs are produced when the electric field reaches the 
critical value E CI along the surface r = r cr . Oppositely charged particles are accelerated 
in opposite directions, leading to the redistribution of the charge bellow the event horizon 
and to a change of the inner geometry of the charged black hole. 



within the event horizon of a charged black hole. The newborn particles do not have 
a spatial infinity to escape to, they are trapped within the event horizon. In the 
Reissner-Nordstr0m manifold in the region between the outer and the inner horizons 
the area coordinate r and the time change their roles. The vector d/dr is now 
timelike, while the vector d/dt is spacelike. The only non-vanishing components 
of the electromagnetic tensor are F uv = —F vu and only F rt = —F tr are nonzero. 
An infalling observer moving along dt/dr = world-line, in the region between the 
horizons, will experience a spatially homogeneous electric field, increasing in strength 
into the future (the r-coordinate decreases). The electric field has only a t-component: 
E = ^P-4:, wherein the bold face denotes usual 3- vectors. The direction of the field 

r z at' 

lines in a regular Reissner-Nordstr0m space-time is from one singularity to the other 
(see Figure |2|). The maximally extended Reissner-Nordstr0m space-time has a charge 
asymmetry in the sense that two external observers in the two past asymptotically 
flat regions, I and I', see the black hole charged oppositely. The left-hand and the 
right-hand singularities seem to such observers to have opposite charges. 

In the interior of a classical charged black hole, between the inner and the outer 
horuzons, pairs of charged particles are produced by the electric field along r = const 
surfaces. An oppositely charged particles are accelerated in the opposite ±d/dt di- 



rections. Thus, if the black hole has a negative charge the positively charged particles 
will be attracted to it, accelerating toward the left-hand singularity, while the nega- 
tively charged particles will be repulsed from this singularity, accelerating toward the 
right-hand one (see Figure ||) . 

This leads to the redistribution of the charge, which was initially concentrated 
near the left-hand singularity. At the end of this process, in the perfectly symmetric 
situation, when the charge is equally divided among the left-hand and the right- 
hand singularities the electric field disappears. But the whole exotic structure of the 
inner part of the analytically extended Reissner-Nordstr0m manifold is due to the 
existence of this electric field. Therefore, vanishing of the electric field leads to the 
disappearance of the tunnels to other asymptotically flat regions. We will show here 
that a similar redistribution of the charge takes place also in a dynamical space-time, 
leading, as we shall see, to the partial closure of the "space tunnels" . 

3 The Physical Model 

In this section we develop an evolutionary formulation which includes the effect of pair 
creation in strong electric fields. We present bellow a simplified toy model, describing 
this effect for a dynamical space-time in which a black hole forms. 

3.1 The Formulation 

The first study of charged particles production in an uniform electric field was un- 
dertaken by J.Schwinger in 1951. The Schwinger's formula for the number of scalar 
pairs created by the field E per unit four-volume is: 

P 2 oo (_-i\n+l ft 

r = d^ D L ^ exp(^L) , (3) 

2iT 2 h c ~^ n z E 

where m and e are the mass and the charge of a created particle and the critical 
electric field, E cr , is defined as: 

eh • [ } 

The production rate of fermions differs from (^ by an overall factor 2 and by the 
absence of the sign interchange (— l) n . 

The effect of the vacuum polarization (the change of the vacuum electric permit- 
tivity) in strong electric fields is instantaneously stronger than the contribution from 



a pair production ( by ln(E/E cr ) for E ^> -E'er), but the latter can accumulate with 
the time |Q. The integrated contribution from the pair creation can dominate the 
contribution from the vacuum polarization. 

Our model utilizes the essential property of the Schwinger's result (|3]) - the ex- 
ponential dependence on the ratio E cr /E, which means that pair creation rate is 
exponentially large for a supercritical field and the rate is exponentially suppressed 
for a subcritical field. This dependence suggests that once the field rises above E cr , 
charged pairs are produced intensively, reducing the field down to the critical value. 

We neglect (i) the net electric current of the pairs and (ii) the stress-energy of 
the produced particles, assuming that the role of the electric current is confined to 
prevent the electric field from rising above the critical value. Hence, the electric field 
will be taken as: 

„ I -'-'ordinary U -'-'ordinary ^ -^cr /p.\ 

^ _C/ Cr II -^ordinary :_ -E'er 

The inordinary stands for the ordinary electric field which would have arised in the 
absence of pair creation. 

In this case we can mimic the effect of particles production as if the system was 
placed in a dielectric medium. Effectively, the polarization of this medium prevents 
the electric field E from rising above E CI , while the electric displacement, D = eE, 
changes. The scalar quantity e is the dielectric constant. The electric displacement is 
related to the density of a free charge via the Maxwell equation: 

V D = 47rp free . (6) 

The dielectric "constant" that leads to (^) is given by: 

1 II -'-'ordinary ^ -C-'cr /_>. 

I *-* | / -C^cr II -'-'ordinary :_ -'-'cr 

The local description of the classical theory of electromagnet ism in a curved space- 
time in a dielectric media can be derived from an effective local Lagrangian: 

£ (cff) = _± E . D = _±rDi!. (8) 

8tt 8tt e V ; 

We use this Lagrangian with e given by (0) to describe the effects of pair creation in 
the strong electric field. 

Our toy model captures the essence of the physics, particularly the fast reduction 
of the supercritical field down to the critical value, E cr and the energy-momentum 



conservation. The dielectric "constant" that we introduce has the same effect as the 
pairs which "shorten" an electric field above E cr . However, we ignore all other features 
of the pairs, specifically we ignore the pairs themselves, their energy momentum tensor 
(which is replaced by a modified electromagnetic energy-momentum tensor that arises 
from the dielectric constant) and their electric current. This is done in order to obtain 
a simple self consistent energy conserving system. It is difficult to estimate what would 
be the effect of the electric current of the pairs. On the other hand it is clear that the 
stress-energy of the pairs would make a positive contribution to the mass-parameter 
(|17|) and will make the effects, which we describe latter, more pronounced. 

Another artificial feature of our model is the ad hoc introduction of E cr . This crit- 
ical field is associated with the mass of the charged particles (see Eq. (^)). However, 
for simplicity our model is based on a massless charged scalar field, whose charac- 
teristics are along null geodesies. For such a massless field the critical electric field 
vanishes and charged massless pairs are produced even for an infinitely small electric 
field. We introduce a critical field E CT as a free parameter which can be used to define 
the mass of the created particles trough the relation (|4]). 

In addition to the above properties there are other minor and physically justified 
assumptions: First, our toy model ignores the contribution from vacuum polarization 
due to the intensive electric field (to be distinguished from the effective polarization, 
which we describe here). This would be justified if the effect of the vacuum polariza- 
tion is small compared to the pair creation contribution. This is in fact the situation 
when the electric field is not exponentially large ||. 

Second, when constructing our model, we have utilized Schwinger's result (|3|), 
more precisely, its exponential dependence on E cr /E. Schwinger's formula is valid, 
strictly speaking, only for an uniform and static electric field over a flat space-time 
background. The approximation of flat-space is valid if the radius of curvature of 
the dynamical geometry is much greater than the Compton wavelength of a created 
particles. The radius of a curvature can be taken of order i? _1//2 , where R is the Ricci 
curvature scalar. The Compton wavelength of the particles is £ — ( ^— J . Thus, 
the condition for a flat-space approximation, £ <C i? -1 / 2 , takes the form 

/ nek \ , , 

R <1. (9) 



But, the Ricci scalar R has been shown, in dynamical models, to diverge approaching 
the Cauchy horizon, (this is the mass inflation scenario) thus, the approximation 
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breaks down in vicinity of the singular Cauchy horizon. Notwithstanding, our model 
is, actually, not based on exact Schwinger result, but on its exponential dependence 
on E cr /E which is non-perturbative. 

Anyway, all approximations will be broken at some moment. Eventually, the 
curvature in the vicinity of the singular Cauchy horizon becomes Planckian since the 
"Coulomb component" of the Weyl curvature diverges exponentially with advanced 
time (for a spherical symmetry- l^l ~ m/r 3 , with m the internal mass parameter.) 
Moreover, the Ricci curvature may dominate the Weyl curvature and surpass the 
Planck values even earlier. In either case our analysis becomes meaningless, and a 
theory of quantum gravity is needed. We obviously consider only the sub-Planckian 
regions. 

3.2 The Equations 

Our goal is to integrate numerically the evolution equations and to follow the collapse 
of a spherically symmetric regular initial scalar field distribution via the formation of 
an apparent horizon and a Cauchy horizon, toward a central singularity. The conven- 
tional choice of coordinates for this dynamical evolution is double-null coordinates. 
In these coordinates: (1) The apparent horizon (when it forms) is regular i.e. it is 
free from unphysical coordinate singularities; (2) For a massless scalar field the char- 
acteristics are null so this choice is "natural" for models involving massless scalar 
fields. 

We choose the line element of the form: 

ds 2 = —a(u,v) dudv + r(u,v) dVL 2 , (10) 

where dQ 2 is the unit two-sphere. There is a coordinate gauge freedom: the choice 
of coordinates u,v is unique only up to a change of variables v' = fi(v),u' = /2(w), 
which leave the line element (10) unchanged. For the time being we do not specify our 



double-null coordinates: these are just general ingoing and outgoing null coordinates. 
Later, when discussing the numerical integration we will fix the gauge freedom and 
specify the coordinates. 

Let the F^ v be the electromagnetic tensor, defined as F^ v = Aw,.^. In a spherically 
symmetric space-time the only non vanishing field components are: F uv = —F vu . 
Thus, only A v and A u need be non vanishing. Then A v can be removed by the gauge 
transformation A a — > A a + A- a , with A = — / A v dv. We are left with A u ^ and we 



denote: a(u, v) = A u . 

We reformulate the set of coupled Einstein-Maxwell-Scalar Field equations as a 
first order system. The numerical integration of the first order system functions very 
well both for the uncharged case, see Ref. [[HJ, and for the charged situation [0. It 



is convenient to define the auxiliary variables: 

d = — , / = r u , g = r v , s = VAiri/;, w = s u , z = s v , (11) 

a 

wherein ip is the complex massless scalar field. We have adopted the notation W x = 

dW/dx for partial derivatives of any function W = W(x,y). 

We denote by q(u,v) the free charge i.e. the charge of the collapsing scalar field, 
and by q the total charge up to the sphere of radius r. The latter is the charge 
defined by the QED effects. The scalar field is collapsing under the influence of the 
total charge q not the free charge q. In a local inertial frame: q = eq and we define: 
a = a/e. 

We write the closed system of equations for the QED-corrected situation: 
Einstein equations: 



El 


= rf v + fg + -a 2 - 


a 2 q 2 

4er 2 


= 0, 




E2 


= 9v~ 


2dg + rz*z = 


= 0, 






E3 


= d u - 


fg a 2 


a 2 q 2 
2er 4 


+ -(wz* 


+ 


Maxwell equations: 














Ml = 


d v - 


a 2 q 
2er 2 


o, 






M2 = 


q v - 


2/ * 


z — 


The scalar field equations: 










Sl 


= rz u + fz - 


f- gw- 


- iedrz - 


f ie 



w*z) H — iea(sz* — s* z) = . 



sz 



16 2 

» H a qs = . 

Aer 

S2 = rw v + gw + fz + iedrz + iedgs + - — a 2 qs = 

Aer 

and, finally, the definitions (pT|): 

D\ ■ 
D2 : 
D3 : 

The QED-corrected equations are reduced to the classical ones by setting e = 1. 
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d- 


Civ 

a 


= 0, 


g- 


r v -- 


= o, 


z — 


s v 


= . 



3.3 The Characteristic Problem 

The system of evolutionary equations have to be completed by a specification of the 
initial and the boundary data along some characteristic hypersurface. For models 
involving massless fields the characteristics are null segments. Thus, it is natural to 
specify the initial conditions and the boundary conditions along null hypersurfaces. A 
satisfactory, for an uncharged case, formulation of the initial-value problem has been 



given by Burko & Ori |16). The generalization for a charged case is given bellow. 

We choose the initial characteristic surfaces to be: the ingoing v = const = v\ 
hypersurface, and, the outgoing u = const = U{ hypersurface. If the domain of 
integration includes the origin of coordinates it leads to the necessity of a series 
expansion of physical quantities in powers of the proper distance from the origin, in a 
vicinity of r = 0. We are, however, interested in the formation of the Cauchy horizon. 
We can, therefore, exclude the origin r = from the domain of integration. We 
achieve this by an appropriate choice of the final outgoing segment u = const = Uf, 
so the domain of integration does not include the origin. 

Now we can remove the coordinates freedom. To do so, we fix the "linear" gauge, 
i.e. we take r to be linear with u or v along the characteristic hypersurfaces. Namely, 
on u = Ui segment we choose g = r v = 1, on v = v t segment we choose / = r u — r u0 . 
To get r along initial surfaces it is necessary to supply r = r(ui,Vi) that serves as a 
free parameter. 

The conventional choice for a characteristic segments Ui = 0, v j = r , yields: 

r(ui, v) = v , r(u, Vi) = ur uQ + r . (12) 

Now we specify freely the scalar field distribution along the initial segments. We 
choose a compact ingoing scalar field pulse along the ingoing u = Ui segment; and 
the "no-perturbation" along the initial outgoing v — v , segment. Specifically, we take 
ipiu, Vi) = that corresponds to a fixed static background for v < Vi. And we choose 
ip(ui,v) = except at some finite region v± < v < i>2, (i>i > v»). To be concrete, for a 
complex scalar field if), ip = <p\ + i<p2, with l5 2 two real scalar fields, we choose: 

V — V\ \ , b . 2 ( V — V\ 



sin In— ,<p2 = —r= sin In , (13) 



An \ v 2 - f i / " v An V v 2 — v x 

where a, b are constant amplitudes, and v' 2 , v 2 are the end-points of each of the real- 
fields pulses, and v± = V{ is their common starting-point. This choice of the initial 
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r-0 singularity 
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of integration 



Figure 3: A schematic diagram of the integrated space-time. 



data is different iable at the matching points v\ 1 v 2 . The integrated space-time is 

schematically depicted on Figure |3|. 

From (|13D we obtain the initial values of z, and w: 

ibir 



z(u i: v) 



W{U,Vi) 



air . / v - Vi 

: sin lii— 



V 2 -Vi 

. 



v 2 -v x 



v 2 -V\ 



sin 2tt 



v — Vi 

v 2 -vi 



(14) 



From the constraint equations E2-E3 and the definition Dl, together with the choice 
a(ui,Vi) = 1 one determines the initial values of d and a: 



d{ui,v) 
+ 



a 2 ii 2 v 



2(v' 2 -v 1 f 

b 2 7T 2 V 



rf2^ + 



sin z 2tt 



\ v 2 -v x 

V — Vl 



v 2 - v x 



(15) 
(16) 



2(v 2 - Vl y 

a{u,Vi) = 1 . 

We assume the Minkowski space-time for v < Vi, therefore, we set q(u, v^ = and 
a(u,Vi) = 0. 

In our coordinates the mass-function (the mass-parameter) becomes: 



.2 4 



m(u,v) = Ul + ^ + —r u r r 



(17) 
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Since the mass-function vanishes for the flat space-time (in the region v < Vj), one can 
calculate: r u0 — — \. It should be noted that for v 3> m our ingoing null coordinate v 
is closely related (proportional) to the ingoing Eddington-Finkelstein null coordinate 
v c . The u coordinate is related to the proper time of an observer at the origin. This 



is defined as [15 



T(u) = / a(u',u')dv! . (18) 

Jo 

In our choice the space-time to the left to the v = v j characteristic hypersurface is 
Minkowskian until the very last moments of the collapse. Hence, a(u,u) = 1 except 
the section when u —* Uf, where a(u,u) —>■ 0. Therefore, the integration in (18]) is 



trivial and yields: T(u) = u. Later, we will find it useful to utilize the proper time of 
an observer at the origin as a measure of the "length" of the Cauchy horizon. 

4 The Numerical Integration Scheme 

We have converted the second order equations to first order equations. Our numerical 
scheme is based on a simultaneous integration of this first order system of coupled 
differential equations. We solve numerically equations El — E2, E4, Ml — M2, SI — 
S2, Dl, D3, .D5. We set e = 1 to obtain the classical collapse with no pair creation. 

The domain of integration is covered by a double-null grid. The characteristic 
initial-value problem is formulated in section |3.3| . The algorithm for the numerical 
integration in the classical case (e = 1) is described in [7J. Here we generalize this 
algorithm to the case when pair creation is included. 

At each step we evolve d and z using E3 and 51 from the hypersurface u to 
u + du. Then we solve the appropriate equations for the rest of quantities along the 
outgoing null rays u + du = const, starting from the initial outgoing hypersurface 
v = Vi. We integrate equation Dl to find a, then we solve the coupled differential 
equations D2 and E2 to get r and g. Next, the equations D3, M2, Ml are integrated 
to obtain s, q and a. Finally, the differential equations El and S2 are solved for / 
and w, respectively. After each step in the w-direction we calculate the electric field 
strength, q/r 2 , along the current outgoing null ray {u + du = const). We use this field 
value to establish the value of e, according to ([7]), for the next it-step. 

This integration scheme uses three distinct methods to evolve the initial data. 
All these methods are well known and commonly used (see, for example, the recipes- 



book by Press W.H., et al. [17]). To evolve the quantities in the u-direction we utilize 
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the 5-th order Cash-Karp Runge-Kutta method. The differential equations in the v 
direction are solved using a 4-th order Runge-Kutta method. The integrations in the 
v direction are performed using a three-point Simpson method. 

It is conventional to define the accuracy of a numerical method by the scaling of 
the numerical error. Thus, nth order accuracy means that, the error scales as the 
step-size to power n: 

F real (x) = F calc (x) + 0(h n ) , (19) 

where F rea i stands for the actual value of a function at a point x, while -F ca i c for a 
calculated at the same point. The Runge-Kutta methods, which we utilize, are all 



at least second order accurate, see Ref. |T7|- The three-point Simpson integration 
method is of order n — 4. 

We have performed a few simulations with the same free parameters, but with 
different grid sizes in order to check the scaling of the numerical errors. The numerical 
scheme proceeds in the u and v directions on grids with corresponding grid-spacings 
h u and h v . These step-sizes are connected by the numerical stability requirement: 
h u < f(u, v)h v , where f(u, v) is a slowly varying function of order unity. We perform 
the convergence test by changing the grids density in both directions. 

In what follows, F h denotes the numerically calculated, with the numerical step- 
size of h, value of a function F at some point. Here h denotes the step-size in the v 
direction. If the numerical scheme, involved in a calculation of F h is convergent, then 
the connection between the above quantities is given by (|19"1). We performed a series of 
numerical simulations with doubled grid densities, or, equivalently, with halved step- 
sizes h v (and, therefore, with halved h u ) . We expect that: F h = F + 0(h n ). Now, 
defining c\ = F h — F h ^ 2 , and c 2 = F h l 2 — F h ^, we expect that: Ci/c 2 = 2 n . Figure 
|](a) depicts log 2 (ci/c2) along an ingoing null v coordinate, for a typical u = const ray. 
It is clear from this Figure that in general n « 3, indicating a third order convergence. 

One, however, notices the point n ~ 0.5, which indicates a very poor convergence. 
Moreover, the plot looks very variable in the region 12 < v < 15. The reason for this 
"jumpy" behavior is understood, if one looks closer at the function F itself. We have 
used F = r(u,v). We sketch on Figure ^(6) the radius r(u,v) along an u = const 
outgoing null ray (see the next section). This sketch is magnified, since the actual 
plots of r for different grid densities are indistinguishable on this scale. In the marked 
box one observes the crossing of curves for different grid densities. This crossing leads 
to a decrease of the convergence order. Notably, after this crossing the convergence 
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(b) 



The lines crossing occur 
in the marked rectangle 
This leads to the loss of 
accuracy. 



(c) 




Figure 4: (a) The error scaling indicates a quadratic or higher convergence except at 
a single point with n — 0.5. (b) To visualize the reason for loss of accuracy we give an 
approximate sketch of r{u, v) as a function of v for u — const, for different grid-densities. 
We give only a sketch and not the real plots, since in the latter the crossing is unobservable 
at these scales, (c) Variation of the calculated r(u, v) as a function of v, for certain value of 
u, for different grid-densities in a typical (with no line-crossing) region. The grid-densities 
are: 60 grid-points per unit interval (solid line), 120 grid-points (dotted line) and 240 
grid-points (dashed line). 
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returns to the high order. 

Figure |](c) displays the variation of the calculated F as a function of v along a 
typical u = const ray for different grids densities. The observed picture confirms the 
convergence. 

5 Results 

5.1 Classical Charged Collapse 

We begin with verifying our numerical code for a classical collapse, a collapse without 
quantum effects. We set e = 1. We fix the free numerical parameters to define the 
problem: a = 0.5,6 = 0.461,^ = r = 10, u/ = 90, v 2 = 16, u^ = 20, e = 0.15. 
The number of grid-points along outgoing and ingoing rays is of order of 10 2 per 
unit interval. We follow the evolution of the regular initial data via the formation of 
an apparent horizon and a Cauchy horizon, toward the central singularity. For this 
specific choice of free parameters the resulting black hole has a charge-to-mass ratio, 
q/m ~ 0.98, in geometrical units. 

Figure |5| displays the metric function r(u, v) as a function of the ingoing coordinate 
v for different values of the outgoing null coordinate u. One can distinguish between 
two kinds of rays: the outer and the inner rays. The outer rays, outgoing rays at small 
u, escape to infinity: r(u,v — ► oo) — > oo. Increasing u, one finds the first outgoing 
ray, which does not escape to infinity, but tends to a constant value. This ray indicates 
the event horizon, which asymptotically coincides with the apparent horizon. The ray 
becomes vertical at about r « 16 (not displayed). Increasing u further, one steps into 
the inner region of the integrated space-time. Outgoing rays in this region approach 
a constant radius, depending on u, in asymptoticly late advanced times v — > oo. This 
indicates the formation of a Cauchy horizon. Following the evolution of the space- 
time, the solution approaches the origin, r = 0. Since our code is not constructed to 
include the origin we stop at u/, before the integration reaches the origin. 

The main feature seen here is the existence of a contracting Cauchy horizon. 
The Cauchy horizon is not a stationary tch — const hypersurface as in the case 
of a classical Reissner-Nordstr0m space-time, but it depends on the outgoing null 
coordinate u, namely, it contracts toward the origin r = in late retarded time u. 

Figure ^ depicts constant radii contour lines in the wf-plane. The bottom of the 
Figure (u = 0) is the initial regular hypersurface. Looking along the u direction 
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Figure 5: Radius r(u, v) as a function of v along u = const rays. The retarded time 
u increases from the initial straight ray toward the origin (from the right to the left), 
indicating the course of evolution of the space-time. 




Figure 6: Lines of constant radius in uu-plane. The radius decreases from the bottom 
toward the top of the Figure. The apparent horizon, along which r v — 0, separates the 
exterior and the interior of the black hole. The latter region contains the Cauchy horizon. 

one soon observes the formation of an apparent horizon, along which r v = 0. The 
apparent horizon separates the two regions with r v > (on the left), and r v < 
(on the right). The latter contains, also, the asymptotically (v — >■ oo) constant, u 
dependent section tqu = const(u), representing the Cauchy horizon. The Cauchy 
horizon itself is a null hypersurface which is approached as v goes to infinity. 
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Figure 7: Logarithm of the mass-function, m(u,v), vs. v along u = const null rays. The 
linear dependence of m(u, v) on late time v indicates the mass-inflation. The retarded time 
u increases toward the top of the Figure. 

We depict the logarithm of the mass-function (|17]) in Figure |7|, for a sequence of 
u = const null rays. The straight lines indicate the exponential dependence of the 
mass-function on v for late advanced times. The exponential growth confirms the 
conclusion that a mass-inflation indeed takes place in this collapse. 

We have performed simulations with different free parameters: changing the am- 
plitudes a and b, the elementary charge e, stretching and squeezing the domain of 
integration. In all these cases we have not seen any qualitative difference between 
the results. The results, which we have obtained, are in excellent agreement with 
the previously established results, and fit well those in M. We conclude, thus, that 
our numerical code gives correct results to this case. We turn now to the problem of 
collapse with pair creation. 

5.2 Collapse with Pair Creation 

The critical electric field, E CT , is an additional free parameter in the problem with pair 
creation. The critical field strength E cr is chosen in the forthcoming graphs so that it 
is reached just after the formation of the apparent horizon. This choice is arbitrary. 
Other comparable values of E CT lead to a qualitatively similar results, unless E CT is 
reached long before the formation of the apparent horizon (see bellow) . 

Figure |3| displays the radius r(u, v) as a function of v for a sequence of u — const 
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Figure 8: The radius, r(u, v), as a function of v along u = const rays for collapse with 
pair creation. The retarded time u increases to the left. One still observes the formation 
of the Cauchy horizon There is also an apparent incline of the late retarded rays toward 
the origin, indicating discharge and signaling the destruction of the Cauchy horizon. 

null rays. This Figure is analogous to Figure [| for the classical situation. On a first 
sight these Figures seem very similar. There are, however, obvious differences, both 
qualitative and quantitative. A closer look uncovers the different incline of the late 
retarded rays (u — > Uf): while these rays depicted on Fig. |5| are practically vertical, 
those on Fig. [8] have an observable incline toward the origin r = 0. The very last 
ray on the latter Figure, has an apparent tendency toward the origin, indicating that 
the r = strong singularity is close. Other signs of the difference between the two 
situations are quantitative ones. The whole "life-time" (in terms of the retarded u- 
time) of the QED-corrected system, before it hits the r = singularity is significantly 
shorter than the "life-time" of the corresponding classical system. For example, for 
the set of free parameters (see at the beginning of the previous section) the evolution 
of the QED-corrected system lasts uj « 12, while for a classical system it takes about 
uj Pd 21, before it crushes into the r = singularity. We note that Uf is the proper 
time of an observer at the origin (|Hf). 

Figure |9| displays the contour lines of constant radii in the wf-plane. This Figure 
is analogous to Figure |^ for the classical situation. Again one can observe a different 
incline of the contour lines in the late retarded and advanced time regions. This 
difference in the incline between the classical and the QED-corrected problems can 
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Figure 9: Contour lines of constant radius in iw-plane for a collapse with pair creation. 
The radius decreases from the bottom toward the top of the Figure. One observes the 
apparent horizon and the Cauchy horizon . 



be interpreted as just a close approach to the intersection of the Cauchy horizon 
with the strong space-like r = singularity during the numerical simulation. The 
straightening of the outgoing rays occurs later in terms of the advanced f-time: the 
curvature is high in the vicinity of the strong singularity turning the rays toward the 
singular origin. The discharge is more apparent in the situation when the black hole, 



as it seen from infinity, has the small charge-to-mass ratio, see Figure 10 



We define the "length" of the Cauchy horizon as a proper time of an observer at 
the origin between the moment when he or she emits a last outgouing light signal 
that escapes to infinity (the event horizon) and the moment when he or she emits a 
first outgoing light signal that unavoidably falls into the spacelike singularity. This 
"length" is related to the retarded time u, see ([TJ). We observe, therefore, that in the 
collapse with discharge the "length" of the Cauchy horizon is "shorter" , compared to 
the classical case. 



On the top panel of Figure 11 we depict the free charge q(u, v) as a function of the 
advanced time v along u = const rays. The bottom panel of Figure |TT] displays the 
total, QED-corrected, charge q(u, v) along the same u = const outgoing surfaces. The 
common characteristic property of the graphs is the straightening of the u = const 
rays in late advanced times, when approaching the formed Cauchy horizon. The 
u = const rays intersect the Cauchy horizon as v — > oo and for different u values the 
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Figure 10: The radius r{u, v) vs. v along a sequence of u = const rays in the situation 
when the charge-to-mass ratio, seen from infinity, of the formed black hole is small compared 
to unity. The bottom panel displays an enlargement of the marked region on the top panel. 
The rays' incline toward the origin is apparent, indicating the discharge. The heavy line 
marks the beginning of the discharge. 



intersection occurs in different points with tch = const{u). Hence, a charge measured 
along an outgoing null rays approaches a constant w-dependent value which decreases 
with u, indicating a contraction of the Cauchy horizon to the r = singularity. 
The difference between the graphs is the strong decay of the charge in the QED- 
corrected case relative to the classical case. In the former one, the charge approaches 
a maximum, which is defined by the strength of the critical field E cr , and then it is 
reduced by the created pairs keeping the electric field at constant E CT value. After the 
pair creation process brought the charge to the "right" value, according to (|7]), the 
charge remains constant along an outgoing u = const rays going toward the Cauchy 
horizon. Moreover, one observes the decrease to nearly zero charge on the bottom 
panel, in contrast to the classical situation. 

The effective dielectric constant e is depicted as a function of the advanced time v 



in Figure [T^. We note the growth of e(u, v) with the lapse of the retarded time u (the 
rightmost ray in Figure is the late u = const depicted ray, the leftmost ray is the early 
u = const ray). This is readily understood: the electric displacement D is growing 
and in order to keep the electric field at the constant E CT value, the polarization field, 
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Figure 11: Top panel: the free charge q(u,v) along u — const null rays. Bottom panel: 
the total charge q(u, v) along the same u = const null rays. On both panels the retarded 
u time increases from the leftmost ray toward the rightmost ray. 

represented by e, must grow, also. 

In Figure [13] we depict the logarithm of the mass- function m(u, v) along a sequence 
of u — const outgoing null rays. The exponential dependence is clear, indicating that 
a mass- inflation takes place also in the collapse with the pair creation. 

In the simulations presented so far E cr was chosen such that E cr was reached just 
after the apparent horizon forms. We have also performed simulations in which the 
critical field was reached (i) long before or (ii) long after the apparent horizon formed. 
Figure |T3](a) displays the space-time of the black hole which is formed in the situation 
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Figure 12: The dielectric constant e(u, v) vs.v along u = const null rays. The retarded 
u time increases from the left to the right. The dielectric constant, which represents the 
"polarization" field, grows with time u. 




135 140 



145 150 155 160 



165 1 70 



Figure 13: The logarithm of the mass- function, m(u, v) vs. v along a sequence of 
u = const null rays. Again, the linear dependence of m(u,v) on the late time v indicates 
the mass-inflation. The retarded time u grows from the bottom to the top. 



(i). In this case the discharge begins long before the black hole forms. The charge 
of the black hole is very small, q <^ m. Therefore r_ « -H-^ -C r + and the Cauchy 
horizon is unobservable, i.e. it forms (if at all) in the domain of very high (Planck) 
curvature. Figure |l4|(b) displays the space-time of a black hole which is formed in the 
situation (ii). In this case the picture is very similar to the classical one. The pair 
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(b) 

Figure 14: The radius r(u, v) vs.v along a sequence of u = const rays in two situations: 
(a) The discharge begins at an early moment of collapse - the Cauchy horizon is unob- 
servable and late rays are falling to the origin; (b) The discharge begins deep inside the 
apparent horizon of the black hole. The Cauchy horizon is almost unaffected relative the 
corresponding classical situation. 

creation begins at a very late stage when a significant part of the Cauchy horizon is 
already formed. The created pairs affect only the last stages of the evolution and the 
full formed Cauchy horizon is only slightly shorter than in the classical situation. 
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6 Summary, Conclusions and Outlook 

We have constructed a dynamical model of a collapse of self-gravitating electrically 
charged massless scalar field including pair creation in strong electric fields. 

Previous studies of the problem were concerned with particles production in the 
electric field: (z) outside 0, [10], [II], |12| or (ii) inside |TJ|, [H|] the event horizon of 
a pre-existing charged black hole. Here, we have formulated the problem in a way 
allowing us to address the question of the influence of the QED-effects on a formation 
of black holes within the framework of an evolutionary model. Our particular interest 
was devoted to a dynamical formation of the Cauchy horizon. 

We have presented and used a toy model that treats the effect of pair creation 
in strong electric fields as an appearance of a local effective dielectric constant. The 
characteristic field strength E cr describing the quantum effects in external field is 
viewed as a boundary between the classical and the quantum stages of the system's 
evolution. We have simulated the collapsing matter by a massless scalar field. Hence, 
this critical field E CT was set as a free parameter which defined the mass of the 
created particles according to @). The conclusions from the numerical integration 
are presented below. 

If the critical electric field strength is reached long before the formation of an ap- 
parent horizon almost complete discharge of the collapsing matter takes place and the 
final black hole is an almost neutral Schwarzschild black hole, see Figure 0(a). The 
black hole's charge, as is measured from infinity is not strictly zero. The remaining 
charge is of order E CT M 2 , where M is the mass of the black hole. 

If the critical value E CT is approached after the apparent horizon forms the black 
hole seems from outside to have all its initial charge. Nevertheless, the process of 
discharge takes place in the inner region of a black hole. In the classical collapse 
strong central Schwarzschild-like singularity, and a weak singular Cauchy horizon is 
formed. With pair creation only a fraction of the Cauchy horizon remains a weak, 
null singular. The rest is replaced by a strong spacelike singularity. The "length" (as 
defined in section (|5.2| ) of this weak singular section of the Cauchy horizon depends 



on the critical field strength E CT . 

It is interesting to consider the interplay between the pair creation effects con- 
sidered here and Planck scale quantum-gravity physics. There are two "extreme" 
situations: (i) The critical field is reached deep inside the inner region, in the Planck 
region surrounding the central singularity. Then the Cauchy horizon is unaffected 
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(or, it is unclear how the horizon is affected at the Planck scales) by the pair creation 
effect. It remains just the same as in a classical collapse (see Fig. |TJ](b)); (ii) The 
critical field is reached at some moment before the formation of the apparent horizon 
then the Cauchy horizon is not formed at all (see Fig. 0(a)). More precisely, as 
it seen by an infalling observer, the formation of the Cauchy horizon occurs on the 
Planck scale. That is, from a practical point of view the Cauchy horizon does not 
form, cf. [IS]. Between two extremes (i) — (ii) the "length" of a weak singular, null 



section of a Cauchy horizon varies between a full, classical "length" and zero. 

It should be emphasized that even if the critical field is approached before the 
formation of an apparent horizon , but after the moment (ii), the weak, null singular 
section of a Cauchy horizon survives. This section is the familiar null singularity, and 
an internal mass-parameter diverges exponentially approaching it (see Fig. [13]). We 
observe in the simulations the formation (or the lack) of a Cauchy horizon in various 
situation, mentioned before. 

A charged spherically symmetric collapse is not a generic phenomena in the nature. 
One does not expect, in general, a significant excess of charge. But even if this 
not a case, the initially charged self-gravitating matter distribution will be rapidly 
neutralized by an accretion of an inter-stellar matter and by the pairs creation process. 

A more generic phenomena of a collapse is one endowed by an angular momentum. 
This is an axisymmetric process and it is rather difficult to study, both analytically 
and numerically. The common property of a charged spherically symmetric collapse 
and the general axisymmetric one is an existence of Cauchy horizons inside the black 
hole. Hence, we use the charged collapse as a toy model to simulate the more general 
rotating situation and to derive a conclusions about the inner structure of the spinning 
black holes. 

It would be interesting question to apply the considerations of this work to a local 
production of particles below the event horizon of a rotating black hole. Outside the 
event horizon of a rotating black hole the coupling of the black hole's spin to the 
orbital angular momentum of particles leads to the phenomena of superradiance. If 
this process takes place below the event horizon of a spinning black hole, then this 
may lead to a destruction of a Cauchy horizon, just like in a charged case. 
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